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£L( Abstract. The Lie point symmetries of the n-th dimensional porous medium 

O ■ diffusion equation are calculated. By restricting to a special class of smooth 

^0 ' functions, the local action of the symmetry group is globalized. This special 

class of functions is constructed using parabolic induction. 



1. Introduction 

The theory of Lie groups finds its genesis when Sophus Lie sets up the task 
to develop an analogue of Galois Theory for differential equations. The original 
prolongation algorithm of Sophus Lie provides a set of point symmetries of a dif- 
ferential equation (for a modern treatment, see [6]). It is well-known that these 
point symmetries generate a local Lie group action on the space of solutions of 
the differential equation and this action seldom globalizes. However, most of the 
modern results on Lie group theory apply to global Lie groups and global Lie group 
representations. As a result, many of the standard techniques of representation 
theory are not always applicable to differential equations. Therefore, it is an im- 
portant problem to find a globalization of the local action of the symmetry group 
of a differential equation. This is usually achieved by restricting to a special class 
of functions. 

In this paper, we study the globalization problem for the n-th dimensional porous 
medium equation, 

u t = A n (u m ), (1.1) 

where A„ is the n-th dimensional Laplacian and u is a function of x = [x\, x n ) G 
K n and tel. When m > 1, this equation models slow diffusion phenomena so this 
condition is often assumed (c.f. [5]). However, for reasons that will become evident 
below, we will allow m S R \ {0, 1}. In fact, we will be particularly interested 
in the case m = *jqr§, because for this special value of m the symmetry group is 
n-dimensions larger than the generic case. The goal of this paper is to describe 
the class of functions on which the action of the symmetry group globalizes and to 
describe the action of the group on these functions. This idea is quite profitable 
when applied to linear PDE's (see [2], [4], and [8]). However, it has also been used 
with non-linear PDE's (see [7]). 

The paper is organized as follows. In Section [2J we list the infinitesimal genera- 
tors of the symmetry group of Equation (jl.l[) . For the one-dimensional case, this 
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is well-known (see pQ or [3]). However, we were unable to find a reference for the 
classification of the symmetry group for n > 2, therefore, calculate it explicitly here. 
Then, we embed the symmetry group as a subgroup of G := SL(2, M.) xSO(n+l, l)o- 
In Section |3l we briefly study the structure of the group G. In Section [4j we use 
parabolic induction to construct a family of representations of G and by restriction, 
of the symmetry group of (jTTTJ) . We use these representations to define class a 
of smooth functions on which the action of the symmetry group globalizes (see 
Equation (|4.ip ). Finally, in Section O we write explicitly the action of the 1- 
parameter subgroups generated by the infinitesimal generators of G. 

2. Symmetry Group 

Using Lie's prolongation algorithm, we calculate the symmetry group of Equation 
(|l.ip ; this yields two cases depending on the value of m. When m ^ the 
infinitesimal generators of the symmetry group are 

Xi = d t (2.1) 
2u 

X 2 = E n + -d u (2.2) 

m — 1 
ii 

X 3 = -td t + -d u (2.3) 

m — 1 

Yi = d t for 1 < % < n (2.4) 

Zi j = x^j — xjdi for 1 < i < j < n (2-5) 

where E n denotes the n-th dimensional Eulcr operator. In this case, they span an 
algebra isomorphic to a parabolic subalgebra of g := s © so(n + 1, 1), where s is 
isomorphic to a parabolic subalgebra of s^IR). 

When m — the symmetry group is generated by the differential operators 
in Equations (12 . 1[) to (I2.5p . together with: 



Wi = (xj - J2 + X ••' V ' + ~~~r d ^ (2-6) 

... ... ^TZ 1 

for 1 < i < n. These infinitesimal generators span an algebra isomorphic to q. The 
isomorphism for the so(n + 1, 1) part is explicitly defined in the following way: 

X2 >-> —E n+ ij l+ 2 — E n+ 2,n+l (2-7) 

^(W t + Yi) h+ E n+lti - E i>n+1 for 1 < i < n (2.8) 

- Yi) ^ E n+2 ,r + E i>n+2 for 1 < i < n (2.9) 

Z i,j ^ E i,i - E 3,i for 1 < i < j < n, (2.10) 

where Ekj is the (n + 2) x (n + 2) matrix with single non-zero entry in the fc-th row 
and Z-th column. There are two possible choices for s; we embed it as the upper 
triangular Borel subalgebra of sb^)- 

3. The Group 

With an eye toward the construction of a family of induced representations of 
the group G := 5L(2,R) x SO(n + 1, l)o, we will study its structure in this section. 
Since most of this material is standard, most details will be omitted (for a general 
treatment see VI of [5]). 
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Let q = {ffip be the Cartan decomposition under the standard Cartan involution. 
Then, t is isomorphic to so(2) x so(n + 1) with so(2) C sbQR) and so(n + 1) 
embedded in so(n + 1, 1) in the upper left (n + 1) x (n + 1) block. The maximal 
abelian subalgebra a C p is given by 
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The centralizcr of a in t, m, is trivial in the st2(R)-component; in the so(n + 1, 1)- 
componcnt, it is isomorphic to so(n) and it is embedded in the upper left n x n 
block. 

livf := E n+lii - E^ n+1 ±E n+2li ±E iin+ 2 e Mat (n+2)x („ +2 ) for 1 < i < n, then, 
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Consequently, the minimal parabolic subalgebras of g are defined asq = mffiaffin 
and q~ = m © a © n~. 

At the level of the group, these subalgebras exponentiate to: 
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The minimal parabolic subgroups corresponding to q and q arc Q = MAN and 
Q~ = MAN~ respectively. 
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4. Induced Representations 

In this section we will construct a family of characters on Q". Then, we will 
induce a representation of G from each of these characters. Via intertwining iso- 
morphisms, we will identify each of these induced representations with a subspace 
of C°°(IR 1 ' r '). From the standard representation theoretical results (c.f. [5]), we 
know that G acts globally on these special classes of functions. Since, the symme- 
try group of the porous medium equation can be embedded in G, the action of G 
descends to a global action of the symmetry group of the porous medium equation. 

We start by noticing that a general element ip G a* is a linear functional that 
acts by ip(H v _ y ) = bv + cy for fixed constants 6,c £ C. By exponentiating ipb -c we 
obtain a character on A. The resulting character, x.r,s '■ A — > C, is determined by 
two continuous parameters r, s S C and defined by 

Xr,s(^a,y) — CL C ^ . 

A character Xp '■ M — > C can be defined by 

Xp(rn JtB ) = (-1)* 

for p£Z 2 . Therefore, we will consider the characters parametrized by two contin- 
uous parameters r, s G C and a discrete parameter p G Z 2 , Xp r s '■ Q ~^ C defined 

by 

Xp,rA<l ) = Xp( m )Xr,s(a) 

where q~ = man~ . 

We consider the induced representation 

Indg_(xp,r,.) := W G C*°°(G) | <p(gq~) = Xp^A^y^id) 

for g G G and q~ G Q - }, 

with the G-action defined by left translation. 

The unipotent group N is isomorphic to M. 1,n via (t, a;) h4 nj iX . Since Q~N em- 
beds in G as an open dense set, it is easy to see from the definition of IndQ_ (Xp,r,s) 
that an element ip G Indg- (Xp,r,s) is completely determined by its value on N. 
Therefore, we can define an isomorphism of vector spaces between Indg- (xp,r.s) 
and the space 

I'(p,r,s) := {/ G G°°(R 1 ' n ) | f(t,x) = <p(n t , x ) for some ip G Indg- (x P ,r,s)} ■ 

(4-1) 

A G- module structure can be given to I'(p,r,s) so that the map tp — > f is inter- 
twining. Therefore, Indg- (x P ,r,s) is isomorphic to I'(p,r,s) as G-modulc. 

5. Actions on I'(p,r,s) 

In this section, we will describe the actions of G and of g on I'(p,r,s). As a 
result, we will find special values of r and s that will determine the class of functions 
on which the action of the symmetry group of Equation (jl.ip . globalizes. 

5.1. Actions of s((2, R) and SL(2, M). In this section, when we consider the action 
of an element g G SL(2, R), it should be understood as the action of its image under 
the natural inclusion map SL(2,M) G. The analogue will be true for elements 
in the Lie algebra. If possible, we will avoid writing this image explicitly to avoid 
cumbersome notation. 
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Proposition 1. An element (" _i) 6 5L(2,K) acts on I'(p,r,s) by 
( a 0a b - 1 ).f(t,x)=sgn(ar\a\ r f(^^,xy 

Let {H,E,F} be the standard basis for s[(2,R). T/ien, i? acts &y f/ie differential 
operator —2td t + rud u and E acts by d t on I'(p, r, s). 

Proof. First, wo write (o-0 = (o\ a )(oa-0 anc ^ we calculate each action of each 
element individually. Let / 6 I'(p, r, s). Then, using the definition of I'(p, r, s) and 
the G action on it, we have 

( o -0 "A*. *) = ( a-0 ■*("*■*) = ^ [( V 2 ) • ^] n t,x)- 

Now we notice that 

[( a Q la a ),In + 2] n t ,x = n o - atiX -[((-- 1 0),J n+2 ]. 
Using the character Xp,r,s and the definition of Indg- (Xp,r,s) ; we obtain 

(S a-0 = X P ,r,,( 

The calculation for the other element is similar, thus omitted. The result on the 
Lie algebra action follows by differentiation. □ 

Remark 1. Setting r = n ^_ 1 we recover the actions of the elements X\ and X$ 
in the symmetry group. We did not calculate the action of F because the space 
of solutions to the porous medium equation is not invariant under the subalgebra 
generated by this element. So, for the goals of this paper, the action of this element 
is not relevant. 

5.2. Actions of so(n+ 1, 1) and SO(n+ 1, l)o- As in the previous section, in this 
section when we consider the action of an element g e SO(n + 1, l)o, it should 
be understood as the action of its image under the natural inclusion map SO(n + 
1, l)o ^ G. The same convention will be used for the actions of the Lie algebra. 

5.2.1. Actions of m and M . 

Lemma 1. Let gij t $ = exp so ^ n+11 ^(9(Eij — Ej^)) for 1 < i < j < n. Then, 

[l2,9i,j,-e] ■ fh,x = n t , gi jt _ ex • [l2,9i,j,-e]- 

Proof. This is a straightforward matrix calculation. □ 

Proposition 2. An element gij.e for 1 < i < j < n acts on a function f G 
I'(p,r,s) by 

9i,j,e-f(t,x) = f(t,g it j-ex). 
The element Eij — Ejs 6 so(n + 1, 1) for 1 < i < j < n acts on I'(p, r, s) by the 
differential operator 

Xi 

Proof. By the definition of I'(p,r,s) and by the action of G on Indg- (Xp,r,s) wc 
have 

9i,j,e-f(t,x) = gi,j,e-f(n t ,x) = f([l2,9i,j,-e]-n t , x )- 
The result now follows from the lemma and the definition of Xp,r,s- The result for 
the Lie algebra follows by taking Jj|e=o- D 



(sgn(o) ( |o j, 1 ^ J , I n +2 ) = sgn(a) p \a\ r f{ a - 
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5.2.2. Actions of n and N. A straightforward matrix calculation shows that 

Tit' ^Xi&i^t^x — nt — t',x — Xiei 

where e; is the standard basis element of R™. We use this fact to prove the following 
proposition. 

Proposition 3. The action ofn t ' iX ' on f 6 I'(p,r,s) is given by 

n t ' <X '.f(t,x) = f(t -t',x- x') 
and vf acts on I'(p, r, s) by the differential operator di, for 1 < i < n. 

5.2.3. Actions of a and A. Using the fact that 

hi,l n t,x = n t ,e'xhi ie , 
we prove the following proposition. 

Proposition 4. The action of h\ lt on f G I'{p,r,s) is given by 

h 1 , e .f(t,x) = e s 'f{t,e*x) 

and i?o,i o-cts on J'(p, r, s) by the differential operator E n + sud u . 

Remark 2. When s = m 2 _ 1 , this action corresponds to the action of the infinites- 
imal generator Xi of the symmetry group of the porous medium equation II) . 

5.2.4. Actions of n~ and N~ . To describe this action, we need to introduce the 
map S : R" — > K dchned by 

S i (x) = l-2x i + \\xf. 
We will also need the map 7i(e, x) : R 1 '™ — ► R™ defined by 

7i(e,x) = (5i(ex) _1 (a; - e||x|| 2 ei) 
and the map Ki : R 1 '" —> 1™ given by 

Ki(e, x) — Si(ex)~ 1 e(ex — e,). 

Then, 

exp(e^r)./(t, x ) = n^ eei .^(n M ) = ^O^o,-^ 

Now, in order to write our result back in terms of / we need to decompose n$ _ ee ..nt >x 
as a product of its N x MAN~ components. 

Lemma 2. For some m 6 M 

Proof. Since SO(n + 1, 1) has real rank 1, W(G, A) has two elements. The non- 
trivial element oj S W(G, A) acts on N by lu~ 1 Nuj — N~. Then, by the uniqueness 
of the Bruhat decomposition, it suffices to show that 

(«t,7 i (e,a;)) n 0,-ee i n *>a:' l l,- log(<5, (ex)) ( n 0,K,(e,K) ) 

differs from the (n + 2) x (n + 2) identity matrix only in the upper left n x n block 
(i.e. it is in M).This is a straightforward but long matrix calculation. 

□ 
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Proposition 5. Let f 6 I' (p, r, s). Then, 

n o,-eeJ(^ x ) = tiifay 'f(t,li(e,x)) 
and v~ acts on I'(p,r,s) as the differential operator 

(xf - x ])di + y)[2 x i x jdj] + 2xsud u . 

Proof. The first equation follows from the previous Lemma and the second equation 
follows form differentiating the first. □ 

We summarize the previous propositions in the following Theorem. 

Theorem 1. The action of the group G on I 1 (p, ^tzj, , f or P £ %2, gives a 

globalization of the action of the local symmetry group of Equation (II. ip . 

Proof. Follows from Propositions [T] to [5] □ 
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